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ABSTRACT 

Magnetars, neutron stars with ultra strong magnetic fields (B ~ 10 14 — 10 15 G), mani- 
fest their exotic nature in the form of soft gamma-ray repeaters and anomalous X-ray 
pulsars. This study estimates the birthrate of magnetars to be ~ 0.22 per century with 
a galactic population comprising of ~17 objects. A population synthesis was carried 
out based on the five anomalous X-ray pulsars detected in the ROSAT All-Sky Survey 
by comparing their number to that of massive OB stars in a well defined volume. Ad- 
ditionally, the group of seven X-ray dim isolated neutron stars detected in the same 
survey were found to have a birthrate of ~ 2 per century with a galactic population 
of - 22,000 objects. 

Key words: pulsar: general — stars: neutron — stars: statistics — stars: magnetic 
fields — supernovae: general 



1 INTRODUCTION 

The end products of massive stars with initial masses in the 
range 8A/q < MinMal < 25M©, neutron stars (NSs) consti- 
tute ~ 0.0001% of the total stellar population of the Milky 



Way in the form of active pulsars ( Lyne & Graham-Smith 
[2006| . It is now understood that not all NSs turn out to ex- 
hibit radio pulsations. This missing link manifests itself in a 



highly energetic class of NSs called magnetars ( Zhang 2003 1 
IDuncan & Thompson ( 1992 1 proposed the existence of mag- 



netars with ultra strong magnetic fields 100 times stronger 
than that observed for normal radio pulsars (B ~ 10 12 G; 



Ziolkowski I2002J) . Magnetar formation, they argued, is 



caused by a vigorous helical dynamo action during the ini- 
tial stages of birth for NSs with initial periods ~ 1 ms. 
Furthermore, the magnetar model predicts the occurrence 
of gamma-ray or X-ray bursts from these powerful sources 
due to starquakes caused by strong magnetic field activity 



( Duncan fc Thompson|1992 Thompson & Duncan 1996). 

Such bursting behavior has been observed for both soft 
gamma-ray repeaters (SGRs) and anomalous X-ray pulsars 
(AXPs; see |Kaspi et al.|2003||Gavriil et al.|2002| ). On March 
5, 1979, a high intensity burst (similar to a short gamma-ray 
burst) and subsequent smaller bursts from the same patch of 
sky recorded by several orbiting spacecrafts unveiled the ex- 



istence of soft gamma-ray repeaters in our galaxy (Mazets 



et al. 19791. Interestingly, the burst from this source only 



lasted ~ 0.1 s, but was followed by a decaying less energetic 
periodic pulse of period ~ 8 s (Terrell et al. 19801. Later 



this source was recognized as SGR 0525-66 in the super- 
nova remnant (SNR) N49, outside our galaxy in the Large 
Magellanic Cloud. Following this discovery, the rest of the 
known SGRs (SGR 1806-20, SGR 1900+14, SGR 1627-41) 
were discovered in our galaxy over a period of 25 years. Sub- 
sequent observations of SGR 1806-20 confirmed that SGRs 
are indeed magnetars characterized by high quiescent X-ray 



luminosities and strong magnetic fields (Kouveliotou et al 



1998) 



The AXPs are another class of very energetic and young 
NSs now believed to be associated to SGRs, and are also 
considered magnetars (sec Mcrcgh etti et al.|2002[ ). As of yet 
only 7 confirmed AXPs have been discovered, namely IE 
2259+586, 4U 0142+61, IE 1048.1-5937, 1RXS J170849.0- 
400910, IE 1841-045, CXOU J010043. 1-721134, and XTE 
J1810-19'f] Both SGRs (in quiescence) and AXPs share 
many similarities and are generally characterized by steady 
X-ray luminosities (L x ~ 10 32 — 10 35 ergs s" 1 ), long spin pe- 
riods (P ~ 5 — 12 s), secular spin down rates with large 
period derivatives (P ~ 10~ n s s _1 ), strong magnetic fields 
(B ~ 10 14 — 10 15 G), and an apparent lack of a binary 
companion ( Woods fc Thompson|2006 1 . Nevertheless, some 
spectral differences exist between the two. Although pulsed 
radio emission has been detected from only one of the mag- 
netars (AXP J1810-197; |Camilo et "aT||2006| |, the rest are 
found to be radio quiet. 

In addition, two new groups of NSs, namely the X-ray 
dim isolated NSs (XDINSs) and the rotating radio transients 
(RRATs), have recently been discovered. The XDINSs are 
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a group of seven radio-quiet NSs (often referred to as the 
'magnificent seven') associated to a nearby OB star complex 



comprising the Gould belt ( Popov et al. 2003 ) . These objects 
have spin periods between 3 - lis, with relatively low X-ray 
luminosities (L ~ 10 30 — 10 31 ergs s _1 ) and magnet ic fields 



(B ~ 10 13 G) compared to the SGRs/AXPs (see Haberl 



2004[|2005||Treves et al.|2000[ ). Here the reader is cautioned 
not to generalize the magnetic field estimate to the whole 
XDINS family as it is solely derived from the spin period 
and period derivative measurements of just two of its mem- 



bers (RX J0720. 4-3125, RBS 1223; |Kaplan fe van Kerkwijk 
2005a|b| . For the remaining XDINSs no period derivatives 
have been obtained so far. As the name suggests, RRATs are 
transient radio sources that exhibit sporadic radio bursts of 
duration between 2-30 ms ( |McLaughlin et al. 12006] ). These 
radio bursts have a maximum flux density of ~ 0.1 — 0.4 Jy 
with burst time intervals between 4 min - 3 hr. What makes 
RRATs worthy of consideration is that their spin periods 
(~ 0.4 — 7 s) and magnetic fields (B ~ 10 13 G) are just in 
the right range to allow suspicion of some association be- 
tween these and the XDINSs ( |Popov||2"o"06| ) . So far 11 such 
sources have been discovered. 

To make the scene even more congested, yet another 
class of NSs begs its introduction. Characterized by long pe- 
riods (~ 0.2 -7s), high magnetic fields (~ 10 13 G), and spin- 
down ages ~ 10 3-4 yrs, high B-field radio pulsars (HBRPs) 
are strong candidates, among the above said objects, in hav- 
ing any connection with the magnetar population ( |Gonzalez| 
et al.|2007 Vranevsevic et aL]|2007 k Many members of this 
class have been detected in the X-rays while others have 
displayed exceedingly dim X-ray emission. 

Due to the paucity of sources, no population synthe- 
sis of magnetars has been performed. As a result, it is still 
unclear how magnetars fit into the grand scheme. Basically, 
the questions that we are asking in this study are: 

(i) What is the birthrate of magnetars? 

(ii) How these objects are associated to other classes of 
NSs? 



Kouveliotou et al. ( 1998 I argued that roughly 10% of all 



NSs become magnetars and that their birthrate is high 
f^Mag ~ 0.1 per century. Estimates of the same magnitude 
have been proposed based on the spin down histories of 
AXPs or derived solely from the age of the SNR to which 
few magnetars are associated ( van Paradijs et al.|1995 |. In 
other instances it has also been noted that spindown ages 
of magnetars may not strictly correspond to the age of the 
SNRs due to decaying fields (jColpi et aL]|2000[|. This ren- 



ders any birthrate calculations to be unsatisfactory since 
these assume constant magnetic fields. Moreover, |Woods fe| 
Thompson ( 2006 1 conjecture that birthrate estimates are de- 



pendent on the efficiency of finding these magnetars which is 
limited by the sensitivity of the detector. Therefore, it is nec- 
essary to conduct a population synthesis in order to obtain a 
reasonable (possibly a lower limit because of transients) es- 
timate of the birthrate of magnetars. Here the authors deem 
necessary to stress the point that this study only looks at 
the persistent AXPs and no transient sources are included in 
the source sample, thus, only a lower limit to the magnetar 
birthrate can be obtained. 

The rest of the paper is organized as follows: Section 
2 outlines the method employed for this study with a brief 



discussion of the ROSAT All-Sky Survey. The results and 
the corresponding discussion are presented in Section 3. Fi- 
nally, we end with a brief note of recommendation in Section 
4 and the conclusion of this study in Section 5. 



2 METHOD 

To provide an overview of the analysis, the whole scheme 
of calculation has been divided into three main tasks. First, 
the detector count rate at a fiducial distance of 1 kpc is 
calculated as a function of the column density Nh, for each 
object in the sample. The sample comprises of the five AXPs 
and the seven XDINSs, all detected in the ROSAT All-Sky 
Survey (RASS) with a high S/N ratio. The SGRs were not 
detected in the statistically well-defined RASS so they are 
not included in the sample. Rather we have simply assumed 
that the existing cadre of instruments could have detected 
any active SGR in the Galaxy. Second, the maximum dis- 
tance to which a source registered a count rate of (a) ^ 0.05 
cts s _1 and (b) ^ Imin (the ROSAT count rate with the 
lowest acceptable S/N ratio) in the RASS, is computed. The 
maximum distance is expected to vary with position on the 
sky due to the difference in the integrated column densities; 
peering into the galactic disk is much harder in comparison 
to the galactic poles. Third, this maximum distance is then 
transformed into a limiting volume within which the num- 
ber of massive OB stars, that are the most likely progenitors 
of neutron stars, is sought. An exercise of this sort requires 
the use of a smoothly varying luminosity function as well as 
the number density of stars in a given volume. Comparison 
of the number of massive OB stars found in the surveyed 
volume to that of their total population yields the scaling 
factor. The scaled populations of AXPs and XDINSs com- 
plemented by the knowledge of their spindown ages, finally, 
provide their birthrates in the Galaxy. 



The ROSAT All-Sky 



(RASS) 



ROSAT, short for ROentgen SATellite, was an X-ray space 
observatory launched on June 1, 1990. Equipped with a 
position-sensitive proportional counter (PSPC), a 2° field 
of view, ROSAT produced an unprecedented all sky survey 
in a period of 6 months. It surveyed the whole sky in the 
~ 0.12 - 2.4 keV (~ 100A - 6A) ene rgy range at a typi - 
cal limiting count rate of 0.015 cts s _1 ( Hiinsch et al.|l999 |. 
ROSAT had an orbital period of ~ 96 min with 40 min night- 
time per orbit and it scanned the whole sky in great circles 
perpendicular to the line connecting the Earth and the Sun. 
Consequently, regions on the sky falling at the ecliptic reg- 
istered much less exposure time (typically ~ 400 s) than 
that at the ecliptic poles (~ 40,000 s). Besides this artifact, 
introduced by the geometry of the survey, each RASS field 
suffered from inconsistent exposure time variations. This is 
illustrated in Fig 1 where the pixels of each field registered 
a wide range of exposure times between 200 s and a little 
over 650 s. 



2.1 Count Rate Vs HI Column Density 

To simulate ROSAT observations of the AXPs in the 0.1 
- 2.4 keV energy range, a two component model (Power- 
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Figure 1. Variable exposure time of an individual ROSAT field. 
The majority of the pixels registered exposure times between 500 
and 650 s. This produced variations in the sensitivity limit of the 
detector for each point on the sky. Thus, the limiting count rate 
is a function of the exposure time. 
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Figure 2. Log of count rate as a function of Nu at a fiducial 
distance of 1 kpc for AXP 4U 0142+61. The sharp drop off is due 
to the exponential attenuation of the source flux with increas- 
ing Nh along the line of sight. Most AXPs, in this study, have 
measured column densities of ~ 10 22 cm -2 . Conversely, XDINSs 
have much smaller column densities due to their nearness. 



law + Blackbody) modulated by interstellar absorption is 
used. Since the XDINSs only exhibit thermal spectra with 
kT < 100 eV ( jPopovpOOo) , only the blackbody component 
is implemented for these objects. 
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(1) 



In the above, r/i is the powerlaw normalization with units 
of photons keV _1 cm _2 s _1 , 772 is the blackbody normaliza- 
tion with units of luminosity per square distance (Lsg/Dio; 
with 1/39 the luminosity in units of 10 39 ergs s _1 and D10 
the distance in units of 10 kpc). With all the spectral pa- 
rameters known the model has to be properly normalized for 
each object. Model normalizations for AXPs 4U 0142+61, 
IE 1 048.1-5937, and 1R XS J 170849. 0-4009 10 are obtained 
from Perna et al. (20011. Normalizations are not available 



for AXPs IE 2259+586 and IE 1841-045, so we use the per- 
centage fluxes of both model components contributing to 



the total unabsorbed flux of the source (|Patel et al. 
2003] ) 



Morii et al. 



2001 



We then normalize the simulated ab- 
sorbed fluxes to that recieved at the telescope to obtain the 
right count rate at the source specific Nh- Since XSPEC 
only provides flux estimates, an ftools subprogram PIMMS 
is used to convert flux into ROSAT count rate. Figure 2 dis- 
plays the log of the simulated count rate as a function of Nh 
at a fiducial distance of 1 kpc for AXP 4U 0142+615 (actual 
distance of source is 3 kpc). In the figure, the exponential 
attenuation of the count rate with increasing column density 
is very much apparent. A sudden change in the slope at low 
column densities is due to the absorption of softer photons. 



2.2 Maximum Probed Distance 

The sensitivity limit of a telescope and the duration for 
which it was exposed to the sky determines how dim, or 
how far, an object could be detected. What is sought here 
is the distance to which a source with known luminosity 
could be detected by ROSAT before its flux drops below 



the detector sensitivity threshold. In this study, two detec- 
tion criteria have been selected. First, a limiting count rate 
of 0.05 cts s -1 (criterion A) is chosen, for it conforms to the 



standards of the RASS bright source catalog (Voges et al 



19961. All five AXPs appear in the RASS-BSC and, thus, 



are well above the standardized detection threshold of a lim- 
iting count rate of 0.05 cts s -1 , a detection likelihoocj^jof at 
least 15, and an accumulation of at least 15 source photons. 
Second, an exposure-time dependent minimum count rate 
(criterion B) is chosen since exposure time variation across 
each RASS field alters the detector threshold in the follow- 
ing way 



1 



+ 



N 



with 



B, 



( 



Bfo\ 



(2) 



(3) 



V Afov - 

In the above, N a = 5 is the signal to noise ratio, B ce u is 
the background count rate in the source detection cell, A ce u 
is the area of the detection cell which, in this study, was 
assumed to be the same as the source extraction area (typi- 
cally ~ 300 arcsec 2 ), Afov is the area of the field of view of 
the telescope, / = 0.9 is the fraction of source counts in the 
detection cell, N m i n = 15 is the minimum number of pho- 
tons required to make a detection, and t is the integrated 
exposure time. The value for Bfov is obtained from the all 
sky background count rate map in the 0.4 - 2.4 keV energy 
range. 

Neutral Hydrogen Model 

The amount of HI integrated along the line of sight limits the 
distance to which a source could be observed. Generally, the 
HI distribution takes the form of an exponentially decaying 
disk, both radially and vertically, away from the galactic 



2 Detection likelihood is defined as 
probability of source existence. 



- ln(l — p) where p is the 
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center. In principle, it should be modified with a warp and an 



additive spiral arm component (Binney & Merrificld 1998). 
For the purpose of this study, a simple exponentially varying 
distribution, of the form given below, is adopted from |Foster| 
& Routledgel (120031). 



n(r, I, b) = n exp 



R(r,l,b) - Ro 



L|b| 



with 

R(r,l,b) 
h z 



R a + (r cos 6) — 2R r cos 6 cos I 
h +h'[R(r,l,b)-R a ] 



(4) 



(5) 
(6) 



where hn = 0.35J? o , h a = 200 pc , b! = 30 pc kpc -1 , and 
-Ro = 8.0 kpc. Notice in (6), the vertical scale height is a lin- 
ear function of the galactocentric radius R. |Foster fe Rout-| 
ledge (20031 argue that for lines of sight probing the very 



midplane of the galaxy, it is proper to adopt such a de- 
pendence as observations confirm a thickening of HI layer 
beyond the solar circle. In order to have a manageable res- 
olution and to keep simulation times shorter, the whole sky 
is divided into 1440x720 points over Right Ascension (RA) 
and Declination (DEC) spanning 360° and 180°, respec- 
tively, with 0.25° intervals. Assuming a reasonable value for 
n = 0.5 cm -3 , the local number density of HI, and inte- 
grating (4) to a heliocentric distance r, for a given galactic 
position (l,b), yields the HI column density. 

N H (r,l,b)= [ n(r',l,b)dr' (7) 
Jo 

The Nh inferred from (7) corresponds to a ROSAT count 
rate from the analysis in the previous section, illustrated 
by Fig 2. Until the count rate fails to satisfy the detection 
criteria, the integral in (7) is iterated with an increasing 
heliocentric distance. Finally, a maximum distance to which 
a given source, in a given direction, could be detected by 
ROSAT is obtained. This exercise is repeated for all twelve 
sources in the sample as well as for both count rate criteria. 



2.3 Starcounts 

In a magnitude limited sample, containing stars brighter 
than a given absolute magnitude M, star counts are usu- 
ally calculated by assuming some stellar luminosity function 
&m(M), and a model of stellar distribution n*(r, I, b) in the 
galaxy. Consider a sample of stars with absolute magnitudes 
between M and M + dM, lying in a specific direction (I, b), 
in a given solid angle ASl, within some distance d. Then 
their total number is simply given by the following expres- 
sion ( |Binney fe Merrifield||1998 1. 

N(Sl,d,l,b)= / $M(M)dM / ru(r,Z,6)r 2 Attdr (8) 
Jm, Jo 



We adopt Mi = —4.9 for the brightest stars and M2 = 
— 1.6 for the dimmest stars that yield neutron stars. This 
particular choice only affects the normalization of the star 
counts but not the derived birth rates which depend on the 
assumed value of n*(r, l,b). For Afi = cosfc Ab Al, where 
Ab = and Al = jfjg, and d = d max from the previous 
exercise, the integral in (8) yields the number of progenitor 
stars in the surveyed volume, towards a given direction. 



Well tested models for the stellar luminosity function 
and the stellar number density are provided in a study done 



by Bahcall & Soneira (19801 (B&S, hereafter). Proposed by 
the authors, their model compares satisfactorily to observa- 
tions made in seventeen different regions on the sky. Thus, 
this study adopts the B&S model with all the underlying 
assumptions. The B&S stellar luminosity function is based 
on starcounts from a subsample of stars with known dis- 
tances and apparent magnitudes. It is derived as a fit to 
observations in the well studied apparent magnitude range 
of +4^m v ^ +22. For -6 < My < 15 



$m(M v 



10 



0{M V -M' v ) 



[l + HT 



( a -f3)S(M v -M')l J 



(9) 



M' v = +1.28, a = 0.74, /3 = 0.04, — = 3.40 



For the number density, the B&S model assumes an ex- 
ponential distribution of stars with all the quantities derived 
from observations in the solar neighborhood. Since stars are 
born out of hydrogen gas, the distribution of massive stars is 
found to very well trace HI in the galaxy, but with different 
scale heights. 



n*(r, I, b) — n+ exp 



R(r,l,b) - Ro 
h R 



Zq + r sin |6| 



(10) 



In the above, n* = 4.03 x 10 pc represents the local 
stellar number density. B&S propose a radial scale length 
for the disk ha ~ 3.5 kpc based on a typical HI model, 
however, this study adopts the same scale length prescribed 
for HI distribution with h R = 0.35i? o = 2.8 kpc. For the 
vertical scale height, B&S rely on the estimates given by 



Schmidt ( 1959 1 which are compounded of many investiga- 



tions, undertaken by various authors (see references therein), 
performed roughly by the middle of the century. Considering 
the availability of new and better data in the past couple of 
decades, it is wise to seek new estimates of the scale heights. 



Such a study was carried out by Reed (2000) where he finds 



a vertical scale height of h z — 45 pc for massive OB stars in 
the solar neighborhood. 

The Gould Belt Enrichment 

A local distribution of massive stars forms a ring like struc- 
ture called the Gould Belt, after B. A. Gould who studied its 
orientation to the galactic disk in 1874. As mentioned previ- 
ously, the seven XDINS which are relatively close by < 400 



pc are associated with the Gould Belt. Grenier (2004 1, in her 



study of the Gould Belt structure, reports that the Sun hap- 
pens to be crossing the belt which is inclined to the galactic 
midplane at an angle of i = 17° with its center 104 pc away 
from the Sun towards the galactic longitude of I = 180°. 
Furthermore, the disk extends to a distance from its center 
of ~ 500 pc and appears to be devoid of any massive stars 
within the inner ~ 150 pc. 

Since the number density of massive stars belonging to 
the Gould belt is roughly twice that of the disk, this enriched 
component must be taken into account while conducting 
starcounts ( Torra et al.|2000 |. Also, it should be noted that 
the Gould belt enrichment is important only for the nearby 
XDINS and does not significantly affect the starcounts in 
the case of the AXPs. This study assumes a constant surface 
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density model for the Gould belt with a Gaussian form for 
the distribution above and below the disk. 



So = 



po exp 



2H 2 , 



dz, 150 pc < d G < 500 pc (11) 



Integrating (11) from — oo to oo yields a relation for the scale 
height 



H., = 



2-Kpo 



(12) 



With the assumption that most stars lie in the disk rather 
than away from it, the surface density can be expressed such 
that S = p h where h is the thickness of the disk. There 
isn't much consensus on what the thickness of the Gould Belt 
is, we adopt a value of ~ 60 pc ( |Grenier|2004[ ). Then, from 
(12) the scale height H z i = 24 pc. The radial distance from 
the center of the Gould disk can be expressed as follows. 



dG{r,l,b) 



(r cos b sin I) 2 + ( 



r cos b cos I + 0.104 



(13) 



+ 



r sin b 



Following the observation by |Torra et al.| ( |2000[ > that p ~ 
2n*, the Gould belt can be incorporated into (10) as such, 



n*(r, I, b) = n* 



x exp 



1 + 2 exp 

R(r,l,b) - R, 
hn 



rsinb x J 
H z > cos i 

Zq + r sin \b\ 



(14) 



for 150 pc < d G < 500 pc. 
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Figure 3. Maximum distance plots for (top) AXP 4U 0142+61 
and (bottom) XDINS RBS 1223. The shading represents the dis- 
tance in units of kpc. Noticeable in the top figure are the disk and 
the central bulge of the Galaxy. This structure is missing in the 
bottom figure due to the low luminosity of XDINSs that makes 
them undetectable even before probing deeply into the galactic 
midplane. The upper limit on the distance in the bottom figure 
lends further support to this fact. In addition, the figures illus- 
trate the sensitivity variations caused by inconsistencies in the 
exposure times as described earlier, and the bright spots mark 
the ecliptic poles. 



2.4 Birthrate 

To deduce the total number of AXPs or XDINSs in the 
Galaxy, the number of sample sources must be scaled ac- 
cordingly. This scaling factor is the ratio of the total number 
of massive stars in the Galaxy to that found in the surveyed 
volume. The former is calculated for a heliocentric distance 
of r — 25 kpc to reach the very edge of the Milky Way. Then, 
the number of massive stars found for each AXP/XDINS are 
averaged and then scaled. 

k 

Naxp = N total ^ (15) 
i=i 1 

where k is the observed population of AXPs. Ni and N total 
represent the number of massive stars in the probed volume 
for each source and the total number of massive stars in 
the whole Galaxy, respectively. Assuming that the popula- 
tions of AXPs and XDINSs are in a steady state, then their 
birthrate is calculated using their spindown ages r c . 



k 

i=i 1 



3 RESULTS AND DISCUSSION 

Maximum distance plots for AXP 4U 0142+61 and XDINS 
RBS 1223 are provided below. Fig 3 illustrates the maximum 
distances to which an object with the spectral characteris- 
tics similar to AXP 0142+61 and XDINS RBS 1223 could 
be detected for criterion B. Strikingly, the curve represents 
the disk of the Galaxy where the gas and dust limits the 
ability of the telescope to peer deep into the midplane. It 
is not surprising to note that the radiation from bright X- 
ray sources is primarily absorbed in the plane of the disk, 
whereas the light travels unimpeded immediately above and 
below it. The knowledge of the maximum surveyed distance 
at each point on the sky enabled us to calculate the num- 
ber of massive, OB stars in the probed volume. For exam- 
ple, starcounts for AXP 4U 0142+61 yielded 267,575 and 
229,967 progenitor stars for limiting count rate A and B, re- 
spectively. Similarly, we found 175 progenitor stars for both 
A and B limiting count rates in the case of the XDINSs. 
Furthermore, integrating the models to the outer regions of 
the Galaxy, we found 523,196 massive OB stars without the 
Gould belt enrichment. Upon its addition, the total number 
increased to 523,694. The final results of this exercise are 
provided in Table 1 with estimates of the population, and 
birthrates, for both classes of NSs. Also listed in the table are 
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Table 1. Listed here for criterion (A) and (B) are the total num- 
ber of objects in the Galaxy, and the corresponding birthrate for 
both AXPs and XDINS. Also, the SGR statistics are provided for 
comparison QKouveliotou et al. [ 19981. 





A 


B 


PaxPs (per century) 


0.20 


0.22 


NaXPs 


12 


14 


PxdinSs (per century) 


2.1 


2.1 


NxDINSs 


22,932 


22,932 


PsGRs (per century) 


~ 0.1 




N SGRs 


3 





the population and birthrate statistics of SGRs calculated 
using (15) and (17) with the assumption of Ni = N to taV, 
however, with so few objects, this value provides only an 
order-of-magnitude estimate. 

Previously, many workers have tried to answer the ques- 
tion asked in this study, and a possible lower limit of 0.1 
magnetars per century has been suggested (Kouveliotou 



et al. 1998 van Paradijs et al. 19951. This study finds the 



birthrate of magnetars to be 0.22 per century, an estimate 
that is in good agreement with the proposed lower limit. 
Essentially, this is the birthrate of active magnetars, in this 
case of active AXPs. If magnetars spend a large fraction 
of the time in quiescence and, consequently, evaded detec- 
tion in the RASS, then the birthrate will be correspondingly 
larger. However, since the birthrate is only a factor of eight 



smaller than the galactic supernova rate (Reed 20051, the 
duty fraction of these objects cannot be too small. In ad- 
dition, from table 1 one will immediately be struck by the 
relatively small difference in birthrates between the AXPs 
and the SGRs. From this we infer like others ( |Zhang||2003| 



Woods fc Thompson|2006 1 that it is possible that the AXPs 
represent a later evolutionary stage of the SGRs. Further- 
more, the birthrate found in this study cannot be assumed 
as the upper limit with utmost certainty due to some of the 
factors described below. 

Firstly, it can be argued that the model assumed for the 
distribution of HI in the galaxy may be too simplistic. It is 
very well known that our galaxy, being spiral in its morphol- 
ogy, has substantial density enhancements along the spiral 
arms. The existence of such a structure can't be ignored, 
yet the model, for this study, has been kept deviod of it. 
Addition of the spiral structure in the HI model would act 
to increase the average sum j- Y2'=i T^F as ^ ne v °l ume wm 
now be less than previously calculated, effectively, reducing 
the number of massive stars in the probed volume. Conse- 
quently, the birthrate estimate will be augmented. Neverthe- 
less, it should be noted that the addition of the spiral arms 
must also increase the number of massive stars as these den- 
sity enhancements are their likely birthplace. After weighing 
both arguments, one can probably conclude that the result 
won't be affected considerably. As an example, a simple error 
analysis for the AXP birthrate, based on the counting error 
of the sample ~ 40% and a naive estimate of the error 

in the spindown ages ~ 50%, amounts to an uncertainty of 
~ 0.2 per century. Similarly, we estimate an uncertainty of 
~ 1 per century for the XDINSs' birthrate. 

Secondly, the form adopted for the distribution of mas- 
sive stars in the Galaxy is similar to that of HI distribution 
with a radial scale height of 0.35-R o . However, the vertical 



scale height was adopted from Reed's ( 2000 1 model. In his 
estimate of the number of massive 03-B2 stars inside the 
solar radius, Reed gives a number of ~ 200,000 stars ( |Reed| 
|2005p . Scaling this number out to the radius of the Galaxy 

with Rg ~ 15 pc by (j^-j , since most of these stars be- 
long to a thin disk, yields a result of ~ 700,000 massive stars. 
Comparing this number to what was obtained for the sim- 
ple model assumed for this study, a conservative estimate of 
523,196 falls short of Reed's estimate. This discrepancy may 
cause a further reduction in the calculated birthrate. 



Thirdly, Colpi et al. (20001 argue that the character- 



istic spindown ages of magnetars inferred from (16) may 
be larger than their true ages. The basis of their argument 
rests on the hypothesis that the magnetic fields, usually as- 
sumed to be constant, are decaying over time. If it is so, then 
a smaller characteristic age would give a larger birthrate. 
Analysis conducted upon fixing for this proposed discrep- 
ancy yields an AXP birthrate of ~ 2 per century. Until the 
issue of the underlying mechanism for magnetar production 
is resolved, it wouldn't be unwise to assume the character- 
istic spindown ages as the true ages of the pulsars. 

This study finds the birthrate of XDINSs to be ~ 2 
per century. Unfortunately, this result was obtained while 
generalizing the characteristic age (~ 1.7 Myr) of just two 
sources (RX J072 0. 4-3125 and RBS 1223; |Kaplan fe van] 
Kerkwijk 2005a|b| to that of the whole sample; furthermore 



Heyl fc Kulkarnrj ( |1998[ ) find cooling ages for these sources 
of about 0.5 Myr. For want of a better determination we 
have taken the mean of these values ~ 1.1 Myr to be the 
typical age of these sources. Under these assumptions, the 
estimated birthrate is in excellent agreement to the value 



proposed by Popov et al. ( 2006 1 for both the XDINSs and 
the RRATSs. 

After having obtained the population size estimates for 
both the AXPs and the XDINSs, it is now possible to find 
any evolutionary affiliation among the different classes of 
NSs, simply by comparison of these statistics. We find that 
there are as many as 2.8 times more AXPs than already 
observed in our galaxy. This corresponds to a population of 
14 AXPs. Including the three SGRs, assuming that all of 
them are observable irrespective of their location, the total 
magnetar population amounts to 17 sources in the Galaxy. 
Similarly, we find that there are ~ 22,000 XDINSs in our 
galaxy. Although the birthrate of XDINSs roughly matches 
the upper limit of that of the normal radio pulsars, the num- 
ber of normal radio pulsars (70,000 - 120,000; |Vranesevic| 
et al.||2004[ ) is larger than that of the XDINSs due to their 



larger typical ages. We find that the birthrate of the AXPs 
and SGRs falls short of the RRATSs (and XDINSs) so that 
they cannot account for the RRAT population unless there 
is a large population of transient AXPs, a similar conclu- 
sion to that found by |Popov et al.| ( |2006[ ) in their analysis 
of magnetar, XDINS and RRAT birthrates. 

On the other hand, the magnetar birthrate is in good 
accord with the HBRP birthrate of ~ 0.2 per century, how- 
ever, the magnetar population again falls short of the HBRP 
population of ~ 187 ± 103 (IVranevsevic et al. 1120071. This 



difference in population simply reuslts from the typically 
larger ages of the HBRPs relative to the AXPs. Consider- 
ing other instances, such as the X-ray emission from many 
HBRPs, the discovery of radio emission from the transient 
AXP XTE J1810-197 ( |Camilo et al.|2006| l, and the compa- 



The Birthrate of Magnetars 7 



rability of spin periods between the two classes, place the 
HBRPs in a circle of suspicion for having some connection 
with the magnetars. 



4 RECOMMENDATIONS 

Recently, hard X-ray emission has been detected from three 
AXPs (4U 0142+61, IE 1841-045, 1RXS J170849. 0-400910) 
by INTEGRAL space telescope operating in the 15 keV - 
10 MeV energy range (Gotz et al. 20061. Furthermore, an 



all sky survey, conducted by using the coded mask IBIS 
telescope on board INTEGRAL, has identified various high 
energy galactic sources, a few among them (~ 42) are yet to 
be identified ( Krivonos et al.|2 007). Considering the results 
of this study and noting the possibility of detecting hard 
X-ray emission from the AXPs, as well as being highly op- 
timistic, we would like to extend our speculation that some 
of the unidentified sources may be those AXPs that have 
evaded detection in the low energy regime. 



5 CONCLUSION 

According to this study, up to 22% of all supernovae explo- 
sions give birth to a magnetar. The small size of the mag- 
netar population reflects their short lifespans during which 
these objects manifest themselves as being highly energetic 
sources in the Galaxy. The birthrates obtained while assum- 
ing decaying magnetic fields, as advocated by |Colpi et al.| 
1 2000), exceed estimates of the supernova rate in the Galaxy 
( Reed|2005 l, thus, this study resorts to constant fields. Fur- 
ther timing analysis is required to gain important insight 
into the characteristics of XDINSs as this has been success- 
fully achieved for only two among the seven. 
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